Late-time tails of a self-gravitating massless scalar field, revisited 
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We discuss the nonlinear origin of the power-law tail in the long-time evolution of a spherically 
symmetric self-gravitating massless scalar field in even-dimensional spacetimes. Using third-order 
perturbation method, we derive explicit expressions for the tail (the decay rate and the amplitude) 
for solutions starting from small initial data and we verify this prediction via numerical integration of 
the Einstein-scalar field equations in four and six dimensions. Our results show that the coincidence 
of decay rates of linear and nonlinear tails in four dimensions (which has misguided some tail hunters 
' in the past) is in a sense accidental and does not hold in higher dimensions. 
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I. INTRODUCTION 



This paper is concerned with the long-time behavior of a spherically symmetric self-gravitating massless scalar field. 
This toy-model of gravitational collapse has been intensively studied in the past leading to valuable insights about 
the validity of the weak cosmic censorship and no-hair conjectures. In particular, Christodoulou proved that there 
are two generic endstates of evolution: Minkowski spacetime for small initial data |JJ] and Schwarzschild black hole for 
large initial data 0] • In both cases the upper bound for the rate of relaxation to the endstate inside the light cone is 
t~ 3 (this was proved in [Ij for the dispersive solutions and recently by Dafermos and Rodnianski Q for the collapsing 
q-i solutions). In view of these rigorous results, one might wonder what is the point of studying this problem again. Our 
I . motivation is twofold. 

^q' First, it is natural to ask whether the decay rates mentioned above are optimal and, if so, what are the corresponding 
i— i, amplitudes of the tails. This kind of quantitative information might be physically relevant (provided that the intuitions 
gained in this toy-model carry over to more realistic situations), for example in assessing the possibility to detect the 
| tails in future gravitational wave experiments. 

• Second and foremost, we want to clarify some longstanding confusion which is widespread throughout the vast 

j relativity literature dealing with wave tails. To explain what this confusion is about we need to make some historical 
remarks. The study of wave tails in general relativity was launched in the seminal paper by Price Q, where he gave 
a heuristic argument that a linear massless scalar field propagating on the fixed Schwarzschild background decays as 

. t near timelike infinity. This result has been later rederived by different methods and confirmed numerically (the 
__ ■ works [1, H, 0, H, [1] are particularly noteworthy), and finally proved rigorously in [|[ (as a special decoupled case 
of the main theorem on the coupled Einstein-scalar system). An especially influential contribution to the study of 

■ tails was made in a pair of papers by Gundlach, Price, and Pullin [g, |lOj. In the first paper, [6[, they argued, using 
linearized theory, that the t~ 3 tail is due to the backscattering of the outgoing radiation off the curvature at large 
• *h , distances, and therefore it is present for any asymptotically flat solution, not only for black hole spacetimes. In the 
second paper, [Io| . GPP solved the spherically symmetric Einstein-scalar field equations numerically and found that, 
indeed, tails do develop for any initial data and moreover they decay as t~ 3 regardless of the endstate of evolution. 
This work was a significant step toward understanding of tails, however the fact that it appeared back to back with 
@ led also (somewhat ironically) to some confusion. Namely, the remarkable agreement between the decay rates of 
tails observed numerically in the nonlinear evolution [To| and the predictions of the linearized theory [f| has been 
interpreted (first rather cautiously by the authors themselves and later with increasing sureness in numerous citations 
of [lCj) as if the linearized theory applies qualitatively (and, as long as the power-law exponents of tails are concerned, 
even quantitatively) in the nonlinear re gime . We wish to point out that this interpretation is too naive. 

We claim that the tails observed in [lfjj| (and later confirmed in [Til Il2| with better numerical accuracy) have 
genuinely nonlinear origin for the dispersive solutions (while for the collapsing solutions they have both linear and 
nonlinear ingredients) . To substantiate our claim we compute the late-time behavior (the decay rate and the amplitude 
of a tail) of the self-gravitating scalar field in even-dimensional spacetime for small initial data using the nonlinear 
perturbative scheme developed in our recent papers [13, LLJ, UM, da, ELD] • The outcome of this simple analytic compu- 
tation is shown to agree extremely well with the results of high-precision numerical integration of the Einstein-scalar 
field equations in four and six dimensions, however it does not agree with the linearized theory in dimensions higher 
than four. Thus, the equality of the decay rates of linear and nonlinear tails seems to be a misleading idiosyncrasy 
of the Einstein-scalar field equations in four dimensions. This paper is concerned only with subcritical initial data 
which lead to dispersion. Work on the collapsing solutions is still in progress. 
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The rest of the paper is organized as follows. In section |TT] we construct a simple iterative scheme for solving 
the spherically symmetric Einstein-scalar held equations in even-dimensional spacetimes. This scheme is applied in 
section IIIII to derive the second-order approximation for the mass function and the third-order approximation for 
the scalar held in four dimensions. The analogous calculation in higher even dimensions is done in section IIVI In 
Section V we compare the nonlinear tails with the linear tails on the hxed Schwarzschild background. Section IVII 
contains numerical evidence confirming the analytic formulae for the tail from sections IIIII and IIVI In section IVIII we 
make some final remarks. 



II. FIELD EQUATIONS AND THE ITERATIVE SCHEME 



We consider the self-gravitating real massless scalar held 4> in d + 1 dimensions, where d > 3 is odd. The Einstein 
equations for the metric g a $ are 

G a f3 = 8nT al3 , T a fj = V^V^ - ^g a0 (V^V^) , (1) 

and the scalar held satishes the wave equation (which is equivalent to V a T Q/3 = 0) 

V M V> = 0. (2) 
We assume spherical symmetry, so <p — <t>(Pi r )i an d use the following ansatz for the metric 

dg 2 = e 2a(t,r) (_ e W{t,r) dt 2 + ^ + ^ (3) 

where dil 2 , _ 1 is the round metric on the unit (d — l)-dimensional sphere. We dehne also the mass function 
m(t,r) = (1 — e~ 2a )r d ~ 2 . In these variables the Einstein equations take the form 

to' = K/ _1 e _2a ^' 2 +e _2 ^ 2 j, (Hamiltonian constraint) (4) 

m = 2nr d ~ 1 e~ 2a cj)(t)' i (momentum constraint) (5) 

= (d-2)-^ T e 2a , (6) 

Stt 

where n — , and primes and dots denote partial derivatives with respect to r and t, respectively. Equation 

d — 1 

(6), corresponding to G\ + G r r = 87r(T t * + TJT) = 0, is sometimes referred to as the polar slicing condition. The wave 
equation which can be viewed as the integrability condition for equations (4) and (5), becomes 

(e-^y--^(r-W)' = 0. (7) 

We assume that initial data are small, smooth, and compactly supported (the last assumption can be replaced by a 
suitable fall-off condition) 

0(O,r)=E/(r), j>(0,r) = eg(r) . (8) 

We make the following perturbative expansion 

m(t,r) = rrio(t,r) + emi(t,r) + e 2 m2(t,r) + . . . , (9) 
(3(t,r) = f] (t,r)+ef3 1 (t,r)+e 2 (3 2 {t,r) + (10) 
<Xi,r) = o (i,r)+e0i(i,r)+e 2 2 (t,r)+£ 3 3 (t,r) + .... (11) 

Substituting this expansion into the held equations and grouping terms with the same power of e we get the iterative 
scheme which can be solved recursively. 

In this paper we consider perturbations about Minkowski spacetime, so too = [3q = cf>o = 0. At the hrst order the 
metric functions toi = Pi = (this follows from regularity at r = 0), while <p\ satishes the hat space radial wave 

equation (□ = d? — d 2 — — — -d r ) 

r 

□0x - , M0, r) = f(r) , <M0, r) = g(r) . (12) 
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The general spherically symmetric solution of equation (JT2J) in odd spatial dimensions d — 21 + 3 is given by a 
superposition of outgoing and ingoing waves 

Mt,r) = ^ et (t,r) + ^ dv (t,r), (13) 

where 

tf-ft r) = — V {2e ~ k)l aik){u) <P*>(t r)- — V( qW W fl4 , 

^ { , ) r z+ l Z^H{l-k)\ (v- u y~ k ' 91 l ' J r£+1 *!(*-*)! (u-u)«- fc ' 1 ' 

and u = t — r, v = t + r are the retarded and advanced times, respectively (the superscript in round brackets denotes 
the fc-th derivative). Note that for compactly supported initial data the generating function a(x) can be chosen to 
have compact support as well (this condition determines a(x) uniquely). 

At the second order Ofa = 0, hence 4>2 = (because it has zero initial data), while the metric functions satisfy the 
following equations 

d~l I 12 , i/2 



m' 2 = KT^ 1 [f 1 +<ff) , (15) 
rh 2 = 2nr d - 1 $ l (j)' 1 , (16) 

« = d») 

We postpone the discussion of this system to the next section and proceed now to the third order, where we have 

□03 = 2/3 2 0i+/3 2 0i+/3^' 1 . (18) 

To solve this equation we use the Duhamel formula for the solution of the inhomogeneous wave equation \3cf> — N(t, r) 
with zero initial data 

t t+r-T 

^ r ) = ^+i/ dr / p e+1 Pi(n)N(T,p)dp, (19) 

|{_ r _ T | 

where Pi(p) are Legendre polynomials of degree i (recall that I = (d — 3) /2) and p = (r 2 +p 2 — (f — r) 2 )/2rp (note that 
— 1 < )Lt < 1 within the integration range). Applying this formula to equation (|18p . using null coordinates rj = r — p 
and £ = t + p, and the abbreviation K(j3, </>) = 2/3</> + $(p + f3'<f)' , we obtain 

i+r t— r 



6>ft r ) = 2*+V+i / ^ / (e-^r +i ^(^(/32(e,?7)^i(e,ry))dr,, (20) 



|t-r| 

,.2 



where now /i = (r 2 + (£— — T]))/ r (^~ V)- If t ne initial data ([8]) vanish outside a ball of radius R, then for t > r + R 
we may drop the advanced part of </>i(t, r) and interchange the order of integration in (|20|l to get 



oo t+r 

1 



fo(*.r)= 2W+ x / dv I ^-vY +1 Pi^)K(M^v)A r i et (^v))dC- (21) 



-oo t—r 



In order to determine the behavior of 4>3(t,r) for large t we need only to know the behavior of the metric function 
@2 (t, r) along the light cone for large values of r (the intersection of the integration range in (|2"Tjl with the support of 
4>i et (t, r ))- This calculation will be done in the next section. Having that, we shall expand the function K in (|21| in 
the inverse powers of (£ — 77) and use the identity (see the appendix in [16j for the derivation) 



:+7 ' dc W = , ^ 2(n - 2Y- r^\t - rjf^ ( 

^ (t ^^r^. ^ ' <^iD 1 1MI U+ ~^2 ~,2ln — 1 



2 ' 2 



v y (2€+ 1)!! [(i - t?) 2 - r 2 ]— 1 I ^ + 3/2 



t-v 



(22) 



where (n - 2)^ = (n - 2)(n - 3) • • • (n - 1 - 1) (^ > 0) and (n - 2)& = 1. 

If one has no fear, this iterative procedure can be continued to higher orders, however it seems like an overkill in 
view of two facts. First, the iteration has no chance to converge (c/.[18|]) so our perturbation series is only asymptotic. 
Second, already the third-order approximation shows excellent agreement with numerical results (see section VI). 
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III. NONLINEAR TAIL IN 3 + 1 DIMENSIONS 



In this section, written mainly for the benefit of the reader who is not interested in higher dimensions, we follow 
the general strategy sketched above to present a detailed calculation of the third-order iterate ^(t, r) in three spatial 
dimensions. In the next section we shall repeat this calculation for any odd spatial dimension d > 3. 

For d = 3 (hence 1 = 0), the solution (|13ll4p of the free wave equation takes the form 



■fc(t,r) 



a(u) — a(v) 



(23) 



where the function a(x) is uniquely determined by the initial data. Substituting (|23[) into (| 1 5j) and integrating, we 
get 



m 2 (t,r) '^4* I [ 2a' 2 (t - p) - d, 



a 2 (t-p) 



dp, 



(24) 



where we used that mi(t,r = 0) = 0, which follows from regularity of initial data at the origin and (|16p . Here and 
in the following we use repeatedly the fact that a (a;) = for |x| > R, R being the radius of a ball on which the 
initial data J3|) are supported. To describe the behavior of m2(t, r) near the lightcone it is convenient to use the null 
coordinate u — t — r instead of t, and rewrite (flM)) as 



/ 7 

ma(«,r) r+ = R Air 2 / a' 2 (s)ds 



a 2 (u) 



(25) 



Next, using the gauge freedom to set /3a(*)?" = 0) = and integrating equation (fl7|) , we get 



fofrr) l = R ^ I 2 



P 

t-p 



a' 2 (s) ds dp - 



a 2 (t - p) 



dp 



(26) 



The first integral can be integrated by parts giving 



r+u>R 



4?r \ ~\\ a ' 2 ( s ) ds + 2 1 



r — (s — u) 



ds 



a 2 {s) 



(r — (s — u)) 3 



ds 



(27) 



In order to determine the tail of 03 we need only two leading terms in the expansion of the above formula in the 
inverse powers of r. Hereafter, it is convenient to define the following integrals (for non-negative integers m, n) 



ds. 



(28) 



Then our results can be cast in the form: 



a , \ r+u>R 47r 

a , \ r+u>R 47T 

P2{u,r) = - — 



2ll(u) + -(2I 2 (u)-I°(u)) + 0^ 
2I° 1 {u) + -{All{u)-a 2 {u)) + o(\ 



r+u>R 47T 



21° (u) - -a 2 (u) + 



(29) 
(30) 
(31) 



Substituting ^ and (|29«3ip into (with t = 0) we get for t > r + 3R 

+oo t-\-r 



&(*,r) = ^ j dr, J 



di 



I°(r 1 )a( V ) +—Q (r l ))+O 



1 



(C-ry) 5 



(32) 



■5 



where 



Qo(v) = 2l\(n)a{n) + (2/ 1 2 (r / ) - 7° ( V ) )a'(n) . 
Elementary integrations over £ and by parts over rj yield for large retarded times u = t — r 



{t 2 - r 2 f 



T + O 



t 2 



where the coefficient 



+00 

T = -2 5 vr J I°{s)a(s)ds 



(33) 



(34) 



(35) 



is the only trace of initial data. From ([34]) we obtain the late-time tails in two asymptotic regimes: <fo(i, r) ~ Tot 3 at 
future timelike infinity (r = const, t — > 00) and (r^>s)(v — 00, u) ~ -r u~ 2 along future null infinity (v = 00, u — > 00). 



IV. NONLINEAR TAIL IN d + 1 DIMENSIONS 

Proceeding along the same lines as described in detail in the previous section, we get the analogues of formulae 
(|29I31[) in d + 1 dimensions (recall that d = 21 + 3 so n = An /(I + 1)) 



a , \ r+u>R {21+ 1)k 
P2(u,r) = 



r 2£+2 



2lj +1 (u) + £ -±l(2lf +1 (u) -(£+ l)7°(u)) + O (i) 

2J0 +1 ( U ) + i±l(4/i +1 ( M ) - (£ + 1) (aW(u)) 2 ) + O (J5) 
2ll 1 (u)-(e+l)i 1 -(a^(u)y+o(± 



■ a 1 s r-+«>fl (2£+ 1)k 
^( u » r ) = ^1+2— 

B'(u r) r+ - >R (2£+1)K 



with I™{u) defined in (f28|) . Substituting these expressions into ([2Tj) we get for t > r + 3R 

+oc t+r 



2 2e + 2 (2£+l) K , 
<h{%r) = — / dn 



1 + 1 



d£ 



-00 £— r 



(f - ?7)^+ 2 



£-77 



+ 7^ i?+i(»7)o W (»7) + T-0/(»7) +0 



d?7 



where 



Q<(»7) = (^ + 2)/| +1 (7 ? )aW(7 ? ) + (27| +1 (^) - (£ + l)7 °(, 7 ))a^+ 1 )(r ? ) 
Using the identity (122|) we get the asymptotic behavior for large retarded times u = t — r 



&(t,r) = 



where 



^2 _ r 2^2^+2 



« + 2-/(l-^)j^ - + | /22 



+00 



= (-l) i+1 2 3i +^ / 4( S )aW( S )d S 



T f 



(36) 
(37) 
(38) 



(39) 



(40) 



-O 



t 



e+i \t 2 -r 2 
(41) 



(42) 



For £ = the formula JUJ) reduces to For £ > 1 the integral can be integrated by parts again, giving 



-f-oc 



t = (-l) i+1 2 3i+5 7r / (a^ +1 \s)) aV-Vtfds. 



(43) 
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Asymptotics at time and null infinity are easily obtained from (|41[) . They read: 

1 



f3£+3 
1 

,,2£+2 



O 



(2€+ 1)!(2£+ 
2(3£ + 2)! 



O 



(44) 
(45) 



This is our main result. We claim that the formulae (44) and (45) provide very good approximations of the tail for 
solutions with sufficiently small initial data. By this we mean that for any given smooth compactly supported profiles 
f(r) and g(r) in (8) one can choose e so small that 



lunt 3e+3 m,r)~e 3 Mt,r)\ = 0(e 5 ) . 

t — *oo 



lim u 2t+2 \{r l+1 (<j)-e 3 fo)){v = oo,u))\ = 0(e 5 ) , 

u — >oo 

at time and null infinity, respectively. Numerical evidence for this claim is given in Section VI. 



(46) 
(47) 



V. LINEAR SCALAR WAVES ON SCHWARZSCHILD BACKGROUND 

For the sake of completeness, in this section we recall briefly what is known about the decay of the massless scalar 
field propagating outside the d + 1 dimensional Schwarzschild black hole 

2 _ ft M \ j + 2 , ft M \ 7 2, „2jo2 



^ = -^l-_j^+^l-_j dr' + r'dnU. (48) 

As above, we consider only odd spatial dimensions d > 3 and use the integer index I = (d — 3)/2. In terms of the 
tortoise coordinate x, defined by dr/dx = 1 — M/r 2l+1 , and the variable ip(x) = r l+1 (f)(r), the radial wave equation 
in the metric (|4"8|) for r > M reduces to the flat spacetimc 1 + 1 dimensional wave equation with the potential (on 
the whole axis — oo < x < oo) 

«f*-<e* + V(x»-0, V-h ")(!^ (49, 



Now, there is an important difference between £ = and £ > cases which is due to the fact that only for £ — the 
tortoise coordinate involves the logarithm. More concretely, for I = we have 

x = r + M\n(r/M - 1) , (50) 

hence for x 3> M 

r = s-MIn(a:/M) + — — - + + 0[ ^ — - , 51 

X X \ x A J 

and therefore 

M ^ 3M 2 \u(x/M) M 2 f M 3 \u 2 (x/M) \ 

which gives rise to the Price tail 0(i,r]_~ Mt~ 3 0]. 
In contrast, for I > 1 we have (see |l5[) 

1 M 21+1 M 2 ( M 3 \ 

r ~ X+ 2l^~ 2£(U + 1) .t«+i + Iv^ 2 " J ' ( } 

which implies that for large x 

tt i , ^+1) (2l+l) 2 (l+l)(4l + 3) M 2 / Af 3 \ 

A remarkable feature of this effective potential is the absence of a term proportional to M. It is for this reason that 
the tail drops very rapidly: 

Af 2 

M»0~^, (55) 

as follows from the general formula t~^ a+2t > for the tail produced by the potential of the form l{t + l)/x 2 + U (x) with 
U{x) ~ x~ a for large x @,[H- 
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VI. NUMERICS 



In order to verify the above analytic predictions we solved numerically the initial value problem ([3][5]) for various 
initial data. To this end we rewrite the wave equation ([7]) as a pair of two first order equations for auxiliary fields 
$ = <// and II = e-' 3 4>: 

if = — L (r d -V$)', and $=(e' 3 n)'. (56) 

We solve these equations with fourth-order accurate Runge-Kutta time stepping using finite differencing in space. At 
each time step we update the metric functions m(t,r) and (3(t,r) by integrating the hamiltonian constraint ((3]) and 
the slicing condition (J6]) with fourth-order Runge-Kutta method. To ensure regularity at the origin we impose the 
boundary conditions $(f, 0) = and II'(i, 0) = 0, which are implemented by 3>(t,r) and U(t, r) being odd and even 
functions of r, respectively. To avoid the contamination of the tail by spurious reflections from the outer boundary 
of the computational grid we place that boundary far away and compute the solution only inside the domain of 
dependence of the initial surface. As was pointed out in Q, a reliable numerical computation of tails requires high- 
order finite difference schemes, since otherwise the ghost potentials generated by discretization errors produce artificial 
tails which might mask the genuine behavior. We used fourth and tenth-order difference schemes for d = 3 and d = 5 
dimensions, respectively 1 . To eliminate high-frequency numerical instabilities we add a small amount of artificial 
dissipation [l9| . that is after each time step advancing solution / from t to t + At on a grid with (At, Ar) mesh we add 
the Kreiss-Oliger dissipative term f(t + At, r) — ► f(t + At, r) + Q k f(t, r), where (for consistency with 2(k — l)-order 
finite difference scheme) Qk is a finite-difference operator of order 2k of the form Qk = (— l) fc+1 (x^) (A+) (A_) 
where a is of order of I and A±f(t, r) = ±(/(t, r ± Ar) — f(t, r)). Finally, to suppress the accumulation of round-off 
errors at late times our codes were run in 128-bit precision. For the above reasons the accurate numerical simulations 
of tails, albeit straightforward, are computationally expensive even in spherical symmetry. 

The numerical results presented here correspond to initial data generated by the function (see (|14p ) 

ea(x) = -J= exp (-x 2 ) (57) 



for different values of e. For these initial data our third-order approximation (44) yields the following asymptotic 
behavior at timelike infinity 

0(t,r) = 12V^£ 3 ^(\ + for d = 3, (58) 

and 

i(,, r) = -?|f!>( 1+0 (!)) for m 

In Fig. I we plot 4>(t, 0) = e^II(t, 0) in d = 3 and d = 5 for three different values of e. The late time tails are clearly 
seen as straight lines on log-log plots. We fit our numerical data with the formula 

r) = At'' 1 exp (B/t + C/t 2 ) , (60) 

which gives the local power index (LPI) [Tl| 

n(t, r) := -t0(i, r)/0(t, r) = 7 + B/t + 2C/t 2 . (61) 

We plot the local power index at r = as a function of \/t in Fig. 2. Our fitting procedure proceeds in two steps. 
First, from the local power index data on the interval < l/t < 1/50 we fit 7, B and C in (|6ip. Next, having 
determined 7, B and C in this way, we fit A in (pIO]) from <fi data on the interval t > 50. We have verified that the 
outcome of the fit (the amplitude A and the decay rate 7) does not depend on the observation point r. The results 
for r — and their confrontation with analytic predictions are summarized in Table 1 for d — 3 and Table 2 for d = 5. 
The agreement between our third-order approximation and the results of numerical integration of the Einstein-scalar 
field equations is excellent for sufficiently small initial data. 



1 On the fixed Schwarzschild background in d = 5 the scalar field <f>(t,r) decays as t~ 10 (see J55}). If a weak self-gravitating scalar field 
decayed at this rate, its tail would be hidden under an artificial tail generated by a ghost potential unless the tenth or higher order 
discretization is used. 





FIG. 2: The local power index n(t, 0) as a function of 1/t. The curves corresponding to small initial data (e = 2 8 , 2 10 , 2 12 ) 
are indistinguishable which indicates that higher order corrections in the perturbation series are negligible. 



In Fig. 3 we plot the fitted amplitude of the tail versus the amplitude of initial data and compare it with our third- 
order analytic formula. The deviation from the scaling A ~ e 3 for large e signals the breakdown of the third-order 
approximation. We stress that we get the same decay rates (7 = 4 for d = 3 and 7 = 7 for d = 5) for all subcritical 
evolutions, regardless of whether our third-order formula predicts correctly the amplitude of the tail (for small data) 
or fails (for large data where higher-order terms in the asymptotic expansion (fTTj) cannot be neglected) . 



TABLE I: d = 3. 



£ 


Numerics: LPI data 


Theory (third order) 


Numerics: <j> data 


B 


C 


7 


7 




A 


2 -iu 


-2.45384 


1.98180 


4.0000 


4 


3.09511e-10 


3.09511e-10 


2 -iu 


-2.44983 


1.71398 


4.0000 


4 


1.98087e-08 


1.98083e-08 


2" 8 


-2.44977 


1.69814 


4.0000 


4 


1.26776c-06 


1.26760e-06 


2-n 


-2.45270 


1.71275 


4.0000 


4 


8.11365e-05 


8.09971e-05 


2 -4 


-2.49938 


1.95863 


4.0000 


4 


5.19274e-03 


5.05355e-03 


2" J 


-2.64286 


2.70811 


4.0000 


4 


0.0415419 


0.0373293 


2~' 2 


-3.14114 


5.24492 


4.0000 


4 


0.332335 


0.222460 


2- 1 


-4.42597 


11.2084 


4.0000 


4 


2.65868 


0.737111 


1 


-6.49635 


18.7893 


3.9999 


4 


21.2694 


1.07316 


2 


-8.98950 


25.4951 


4.0002 


4 


170.156 


0.997247 


4 


-11.8828 


32.3576 


4.0021 


4 


1361.24 


0.892345 
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FIG. 3: The log-log plot of the amplitude of the tail as a function of the amplitude of initial data (black dots). The third-order 
approximation (dashed line) is excellent for small data, but it deteriorates for large data lying near the threshold for black hole 
formation (e ~ 1). The scaling A ~ e 3 was previously observed by GPP (see Fig. 14 in [Ty| ) . 



TABLE II: d = 5. 



e 


Numerics: LPI data 


Theory (third order) 


Numerics: cj> data 


B 


C 


7 


7 


A 


A 


2 _ia 


4.64290e-04 


0.650257 


7.0000 


7 


-1.21667e-08 


-1.21667e-08 


2 -uj 


2.48554e-04 


0.654252 


7.0000 


7 


-7.78672e-07 


-7.78644e-07 


2 -8 


-1.02440e-04 


0.654138 


7.0000 


7 


-4.98350e-05 


-4.98072e-05 


2 -« 


-6.10626e-03 


0.660806 


7.0000 


7 


-3.18944e-03 


-3.16123e-03 


2-4 


-0.100262 


0.765030 


7.0000 


7 


-0.204124 


-0.177845 


2" 3 


-0.380511 


1.07174 


7.0000 


7 


-1.63299 


-0.986771 


2~' 2 


-1.28878 


2.05237 


7.0000 


7 


-13.0639 


-2.82527 


2" 1 


-3.48234 


4.29513 


7.0000 


7 


-104.512 


-3.53701 


1 


-6.86634 


7.06840 


7.0000 


7 


-836.092 


-3.25661 



VII. FINAL REMARKS 



Using the third-order perturbation method we derived explicit formulae for the late-time tail (the decay rate and 
the amplitude) of a spherically symmetric, self-gravitating massless scalar field for solutions starting from small initial 
data. We verified that these formulae are in excellent agreement with the results of numerical integration of the 
Einstein-scalar field equations in four and six dimensions. Our results show that the tail has genuinely nonlinear 
origin and should not be mistaken with the linear tail coming from the backscattering off the Schwarzschild potential. 
It seems to us that this distinction between linear and nonlinear tails has not been widely recognized in the past 
which is probably due to the fact that in four-dimensional spacetimes these two different tails decay at the same rate 
t~ 3 . To demonstrate that this coincidence is an idiosyncrasy of four dimensions, we computed both kinds of tails in 
d + 1 dimensions for d = 2£ + 3 (£ = 1, 2, ...) and showed that the linear and nonlinear tails decay at different rates: 
£-(6£+4) anc j £-(3£+3)^ respectively. This illustrates how viewing the dimension of a spacetime as a parameter may 
help understand which features of general relativity depend crucially on our world being four dimensional and which 
ones are general. 

It would be interesting to generalize the results of this paper to collapsing solutions where the endstate of evolution 
is a black hole. The studies in this direction are in progress and will be reported elsewhere. We expect that in 
this case the tail has both linear and nonlinear contributions with the latter being qualitatively the same as for 
dispersive solutions described above. Note that the analogous perturbative calculation of tails is much harder on the 
black hole background because two basic tools that we used above, Huygens' principle and the explicit expression 
for Duhamel's formula, are missing. For dispersive solutions these tools allowed us to compute the third-order 
perturbation in a shamelessly explicit way, however, from the perspective of generalizing the results to collapsing 
solutions, it is instructive to redo this calculation in an asymptotic manner keeping track of only leading order terms 
in the perturbative equations. Such an efficient calculation, which gives additional insight into the mechanism of some 
cancelations in our asymptotic expansions, has been done recently by Szpak [20| and will appear as a comment to 
this paper. 
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